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ABSTRACT 

Conservation of the phase-space density of photons plus Lorentz invariance requires 
that the cosmological luminosity distance be larger than the angular diameter distance 
by a factor of (1 + z) 2 , where z is the redshift. Because this is a fundamental symme- 
try, this prediction — known sometimes as the "Etherington relation" or the "Tolman 
test" — is independent of the world model, or even the assumptions of homogeneity and 
isotropy. It depends, however, on Lorentz invariance and transparency. Transparency 
can be affected by intergalactic dust or interactions between photons and the dark sec- 
tor. Baryon acoustic feature (BAF) and type la supernovae (SNela) measures of the 
expansion history are differently sensitive to the angular diameter and luminosity dis- 
tances and can therefore be used in conjunction to limit cosmic transparency. At the 
present day, the comparison only limits the change At in the optical depth from redshift 
0.20 to 0.35 at visible wavelengths to At < 0.13 at 95% confidence. In a model with 
a constant comoving number density n of scatterers of constant proper cross-section a, 
this limit implies na < 2 x 10~ 4 h Mpc -1 . These limits depend weakly on cosmological 
world model. Assuming a concordance world model, the best-fit value of At to current 
data is negative at the 2a level. This could signal interesting new physics or could be the 
result of unidentified systematics in the BAF/SNela measurements. Within the next 
few years, the limits on transparency could extend to redshifts z ~ 2.5 and improve 
to < 1.1 x 10~ 5 h Mpc -1 . Cosmic variance will eventually limit the sensitivity of 
any test using the BAF at the na ~ 4 x 10~ 7 h Mpc -1 level. Comparison with other 
measures of the transparency is provided; no other measure in the visible is as free of 
astrophysical assumptions. 

Subject headings: cosmology: observations — cosmology: fundamental parameters 
- large-scale structure of universe — radiative transfer — relativity — supernovae: 
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Introduction 



The transparency of the universe is extremely good. A typical astronomical camera has a 
shutter whose thickness is measured in microns; that shutter is far more opaque than the entire line 
of sight to the majority of extragalactic sources, even at extremely high redshifts, despite — in many 
cases — considerable column depths of dark matter, plasma, gas, and dust. There are, however, very 
few quantitative measures of the transparency with contemporary astronomical data. 

There are several sources for photon attenuation that are clustered with matter. For example, 
as stars eject heavy elements, they also eject photon-absorbing ash (called "dust"). The gas and 
plasma in and around galaxies absorbs, scatters, and re-emits at longer wavelengths some fraction 
of incident radiation. More speculatively, if the dark matter is an axion or axionlike particle, it will 
in general have photon interactions, which can in principle produ ce effective absorption of photons 



in re gions of high dark matter density and high magnetic fields (jSikivie 



1983 



Raffelt fc Stodolskv 



19881 ) . The sources of attenuation — such as these — that are clustered with matter will be correlated 
with galaxies and large-scale structure, and can be found with "angular difference" measurements 
that compare the apparent properties of sources whose lines of sight have different impact param- 
eters with the correlated structure. 

The Sloan Digital Sky Survey has permitted very sensitive angular difference measurements, 
which find that the attenuation correlated with large-scale structure is very small and consistent 
with being caused by dust, presumably the dust emitted by the stars in the galaxies that populate 
the struct ure. Measurements in the literature constrain this in visible bandpasses at the part in 
10 3 level (jMenard et aZ.ll2Q08l : LBovy et al.ll2008l ). To be specific, these studies constrain differences 
in opacity along different lines of sight caused by absorbers correlated with galaxies. 

It is possible, however, that there might be unclustered or "monopole" sources of attenuation, 
that affect all lines of sight equally, for example if the nonmatter contributors to the cosmological 
energy— momentum tensor (the "dark energy" in modern parlance) have interactions with photons, 
or if there are small violations of Lorentz invariance on cosmological scales. These sources of 
attenuation are much harder to detect with differential experiments, but they can be detected by 
"radial difference" experiments that compare cosmological sources of radiation of known physical 
properties at different redshifts or radial distances. 

A number of different mechanisms have bee n proposed during the last decade to ex plain the 
observed dimming of type la supernovae (SNela; iRiess et aZ.lll998l : iPerlmutter et a/.lll999i ) without 
cosmic acceleration by employing exactly such unclustered sources of attenuation. The mixing of 
photons with axions in extragalactic magnetic fields could lead to photons oscillating into axions 
with a non-negli gable probability over cosmolog ical distances, thus reducing the flux of SNela at 
large distances (ICsaki. Kaloper &; Terningl 12002). Alternatively, "gray" intergalactic dust could 
be so gray as to evade detect ion through it s reddening, while still being cosmologically important 
because of its overall opacity (|Aguirreill999l ). In order to account for the observed SNela dimming, 
these models predict violations of transparency at the order-unity level out to redshifts of unity 
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(e.g., iMortsell. Bergstrom. Goobarll2002l ) 



Furthermore, even if there are no exotic absorbers in the universe, it is difficult (and usually 
model dependent) to infer the total mean opacity from any absorbers that have been found by 
angular difference experiments. Radial difference and angular difference experiments are comple- 
mentary, even when the absorbers are mundane; although radial difference experiments are usually 
less precise, they provide irreplaceable information for measuring total opacity. 

Radial difference ex periment s are so metimes known as "Tolman tests" because they are variants 
of the test proposed by Tolman ( 193ol ) of the expansion of the universe: a test of the (1 + z)~ 4 
(where z is the redshift) dimming of bolometric intensity (energy per unit time per unit area per 
unit solid angle; also called "bolometric surface brightness") with redshift. The intensity is closely 
related to the phase-space density of photons, which is conserved (in a transparent medium) along 
the light path; that conservation plus Lorentz invariance implies the (1 + z)~ 4 relation: one factor 
of (1 + z)~ l comes from the decrease in energy of each photon due to the redshift, another factor 
comes from the decrease in photons per unit time, and two more factors arise from the solid-angle 
effects of relativistic aberration. The Tolman test does not really test for the expansion of the 
universe — the result does not depend on cosmological model, or even assumptions of isotropy or 
homogeneity — but rather for the combination of conservation of photon phase-space density and 
Lorentz invariance. 

In addition to models that violate transparency, there are models that violate Lorentz in- 
variance. Generically these models produce an energ y-dependent speed of light and birefringence, 



break ing the perfect nondispersiveness of the vacuum (|Amelino-Camelia et al 



1998; Gambini fc Pullin 



1999). These effects generally become larger with increasing energy, and observations of high-energy 



sources such as active galactic nuclei (IBiller et al. 



and gamma-ray bursts (GRBs; ISchaefer et al 



1999 



1999 



Ellis et al. 



Aharonian et a/.M20Q8l : lAlbert et a/.l l2008) 



2006) have shown that the linear 



dispersion relation for photons is preserved to good accuracy at these energies. Therefore, while 
these models do fail the Tolman test because of their nontrivial dispersion relations, the effect will 
be unmeasurably small for low-energy (visible-band) photons. 

By far the most precise radial difference test to date has been performed in the radio with the 
cosmic background radiation. In contemporary cosmological models, the CBR comes from redshift 
~ 1100 and is a near-perfect blackbody. The COBE DMR experiment established that the spectrum 
and amplitude of t his radiation are co nsistent with the blackbody expectation at the < 10~ 2 level 
at 95% confidence (IMather et aZ.lll994l ) . A source of attenuating material, unless in perfect thermal 
equilibrium with the CBR, would tend to change either the spectrum or the ampl itude, so this 



resul t provides a very strong constraint on the transparency at cm wavelengths (IMirizzi et al 



2005)). Another test of transparency at cm wavele ngths is the inc r ease in the CMB temperature 
?cmb according to the relation Tcmb oc (1 + z). ISrianand et al\ (|2000l ) find consistency with a 
transparent universe by measuring the CMB temperature at z = 2.3. Of course, many sources of 
attenuation are expected to be wavelength dependent, so these beautiful results may not strongly 
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constrain the opacity in the visible. 



At visible wavelengths, there have been much less precise radial difference tests performed with 
galaxies, whose properties would deviate from naive predictions under extreme attenuation. After 
correcting for the evolution of stellar populations in galaxies, these studies find consistency wit h 
transparency at the 0.5 mag level at 95% confidence (jPahre et q/Jll99a ; iLubin &: 
which correspond to optical depth limits < 0.5 out to redshift z ~ 1. Unfortunately, the precision 
of these tests is not limited by the precision of the measurements, but rather by the precision 
with which the evolution of galaxy stellar populations is known; the results will not be improved 
substantially with additional or more precise observations. 

Another test of transparency at visible wavelengths involves the measurement of the Cosmic 
Infrared Background (hereafter CIB). The absorption of visible photons by a diffuse component of 
intergalactic dust and its re-emission in the infrared contributes to the CIB. The amount of dust 
required to explain the s ystematic dimming of SNela would produce most of the observed CIB 
(jAguirre &: Haimanll2000l ). However, discrete sources (e.g., dusty star- forming galaxies) also emit 
in the infrared and account for almost all of the CIB, strongly constraining the role of dust in the 
dimming of SNela. An y constraint on the tran sparency from the CIB requires a careful subtraction 
of the discrete sources (IHauser h Dwekll200ll ). 



The Tolman test can be rewritten as a relationship among cosm ological dist ance measures. 
There are several empirical definitions of distance in cosmology (e.g., iHoga Il999i ); the most im- 
portant for contemporary observables are the luminosity distance Z?l and the angular diameter 
distance Da- The luminosity distance -Dl to an object is defined to be the distance that relates 
bolometric energy per unit time per unit area S (flux) received at a telescope to the energy per 
unit time L (luminosity or power) of the source, or 

L 



4tt Dl 



(1) 



The angular diameter distance Z?a is the distance that relates the observed (small) angular size 
measured by a telescope to the proper size R of an object, or 



G = — 



(2) 



Because the ratio of flux to the solid angle is essentially the intensity, the (1 + z) 4 redshift 
dependence of the intensity is reflected in these distance measures by 



D L = (1 + zY D A 



(3) 



Both distance measures are strong functions of the world model, but this rel ationship between 
them — known sometimes as the "Etherington relation" (after lEtheringtonlll933l . who showed that 
the result is valid in arbitrary spacetimes) — depends only on conservation of phase-space density of 
photons (transparency) and Lorentz invariance. Fortunately, for some fortuitous types of objects, 
these distances can be measured nearly independently. 
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A test of this type for transparency has been proposed and carried out previously (jBassett fc Kunz 

2004al fbh. with luminosity distances from SNela and angular di ameter distances estimated from 



FRIIb radio galaxies, compact radio sources, and X-ray clusters (jUzan et aZ.ll20Q4l ; I Jackson! 120081 ) . 
The results were imprecise because there are many astrophysical uncertainties in the proper diam- 
eter estimates of these exceedingly complex astrophysical sources. 

In the contemporary adiabatic cosmological standard model, there is a feature in the dark- 
matter autocorrelation function (or the power spectrum) corresponding to the comm unication of 
density perturbations by aco ustic modes during the period in which radiation dominates ([Peebles &: Yu 
1970; Eisenstein et aZ.1 120051 ) . This feature has a low amplitude in present-day structure (that is, 
the distribution of galaxies), but because it evolves little in comoving coordinates, it provides a 
"standard ruler" for direct measurement of the expansion history. A measurement of the baryon 
acoustic feature (BAF) in a population of galaxies at a particular redshift provides a combined 
measure of the angular diameter distance to that redshift (from the transverse size of the fea- 
ture) and the Hubble constant or expansion rate at that redshift (from the line-of-sight size of 
the feature). As we discuss below, as signal-to- noise improves, the BAF can be used to measure 
the angular diameter independently of the local Hubble rate. Most importantly, because the BAF 
arises from extremely simple physics in the early universe when the growth of structure is linear and 
electromagnetic interactions dominate, the BAF measures the angular diameter distance with far 
fewer assumptions than any method based on complex astrophysical sources in the highly nonlinear 
regime. 

At the same time, SNela have been found to be standard — or really "standardizable" — c andles. 



1938; 


Tammann 


1979; 


Colgate 


1979; 


Riess et al. 


1998; 



Perlmutter et aZ.I Il999l ) . Up to an overall 



scale and some uncertainties about the intrinsic spectra and variability among SNela, a collection 
of SNela measure the luminosity distance. 

Given overall scale uncertainties, the most robust test of global cosmic transparency that can 
be constructed from these two distance indicators is a measurement of the ratio of the distances to 
two redshifts z\ and z%. That is, transparency requires 



D h (z 2 ) [1+Z 2 ] 2 D A (z 2 ) 



D L ( Zl ) [1 + zi] 2 D A ( Zl ) 



(4) 



This expression cancels out overall scale issues and is independent of the world model. We perform 
a very conservative variant of this test below, where we measure the left-hand side with SNela and 
the right-hand side with the BAF, marginalizing over a broad range of world models. 

The tests presented here are not precise, simply because at the present day BAF measurements 
are in their infancy, and we make use of no cosmological data other than the BAF and SNela. As we 
discuss below, when these measurements are made at higher redshifts and with higher precisions, our 
limits on transparency and Lorentz invariance will improve by orders of magnitude. Eventually they 
may be limited not by the data quality but by the cosmic variance limit on the BAF measurement 
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itself (|Seo Eisensteinl 120071 ). 



2. Data, procedure, and results 

In surveys to date, where the BAF is measured at low signal to noise, the optimal extraction of 
the signal best constrains not the angular diameter distance directly, but rather a hybrid distance 



D 



v 



cz[l + z] 2 D\ 
H{z) 



1/3 



(5) 



where Z?a is the angular di ameter distance and H(z) is the Hubble expansion rate (velocity per 
unit distance) at redshift z (jEisenstein et aZ.ll2005l ) . 



Using data from the Sloan Digital Sky Survey and the Two Degree Field Galaxy Redshift 
Survey, the power spectrum and BAF have now been measured in samples of massive, red galax- 
ies at two different redshifts: z = 0.20 and z = 0.35. The measured BAF at each redshift z 
translates to a distance measure Dy(z). Accounting for covariances in the measurements at the 
two redshifts (which are not based on entirely inde pendent data sets) , the ratio of distances is 
Z?v(0.35)/L> v (0.20) = 1.812 ±0.060 (68% confidence: [percival et alhooj ). 



We formed two samples of SNela data from a recent compilation (jDavis et all 120071 ). "Sam- 
ple A" consists of all seven SNela in the redshift range 0.15 < z < 0.25 and "Sample B" consists 
of all 22 SNela in the redshift range 0.3 < z < 0.4. We estimate the distance— modulus DM at 
z = 0.20 and z = 0.35 by fitting a straight line to Samples A and B separately (Figure [JJ, and 
obtain a distance— modulus difference 



AD M bs = DM ohs (0.35) - DM obs (0.20) = [1.34 ± 0.09] mag 



(6) 



where we are indicating that this is an observed value, and might differ from the true value if there 
is opacity. 

The distance modulus derived from the SNela is systematically affected by the presence of any 
intervening absorber. Let r(z) denote the opacity between an observer at z = and a source at 
redshift z due to such extinction effects. The flux received from this source is reduced by the factor 
e — t(z)_ 'pkg m f errec i (observed) luminosity distance differs from the "true" luminosity distance: 

D L 2 ohs (z) = D L i nc (z)e^ . (7) 

The ratio of the luminosity distances at two different redshifts Z\ and z-i depends upon the fac- 
tor e [ r ( - 35 )- r ( - 20 )]/ 2 . xhe inferred (observed) distance modulus differs from the "true" distance 
modulus: 

DM ohB (z) = DM tvue (z) + [2.5 loge] r(z) . (8) 
Taking differences of distance moduli at the two redshifts: 



ADM ohs = ADM true + [2.5 loge] At 



(9) 
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where At = \t{z2) — r(zi)]. If the distance indicator from the BAF is unaffected by the absorption 
as we expect, then 



At 



ln(10) 
2.5 



ADMob , _ 7 . 5 log + ,5 ,og ( g|±i&j 

\Dy(zi)J \z 1 [l + z 2 \ 2 H(z 2 , 



(10) 



The above equation can be used to determine At from z = 0.35 to z = 0.20 in light of the 
ratio of the distances Dy obtained from the BAF observations (hereafter B) and the difference in 
distance moduli obtained from the SNela observations (hereafter S) at these redshifts. However, 
the last term in the above equation makes the result cosmology dependent. Therefore, we follow a 
Bayesian approach and assign posterior probabilities to 100 uniformly spaced values of At £ [0,0.5] 
by marginalizing over 100 x 100 ACDM cosmologies uniformly spaced in the (£I a ,£Im) plane with 
ft A G [0,1] and Q M G [0,1]. Thus, 

p{At\s,b)= [ f p(n A ,n M \B)p(A T ,n A ,n M \s)dn M dn A , (11) 

where P(Q a ,Qm\B) and P(At, Q a , Qm\S) are the posterior probabilities of the set of model pa- 
rameters given B and S respectively. We assume that the uncertainties on B and S are Gaussian 
and calculate the likelihood of B and S for different sets of parameters in the (At, Qa, ^m) space. 
Assuming flat priors on ft a and ft a/ in the ranges < ft < 1, and flat prior on At in the range 
< At < 0.5, the posterior probabilities P(Q a ,Qm\B) and P(At,Q a ,£Im\S) are calculated from 
the likelihoods of the two data sets. Equation (jlll) yields the posterior for At, marginalized over 
all world models. Figure [2] shows the posterior P(At|5, B) for the difference in optical depths 
between redshifts 0.35 and 0.20 obtained from the procedure outlined above. The posterior peaks 
at and yields At < 0.13 at 95% confidence. The result demonstrates the transparency of the 
universe between these two redshifts, although not at high precision. 

The abundance and absorption properties of absorbers can be constrained using the difference 
in optical depths measured above. Let n(z) denote the comoving number density of absorbers, each 
with a proper cross-section 17(2;) at redshift z. The difference in optical depths between redshifts 
Z\ and Z2 is then given by 



E(z) 



where Dn is c/Hq and 



f Z2 (l+z) 2 
At= n(z)a(z)Dn 1 ZJ dz » ( 12 ) 

J Zl 



H(z) 



E(z) = — ^ = vAMl + zf + n fe (i + z) 2 + fi A • (13) 

Bq 

In detail, the output of this integral depends on world model. For the concordance model, Hubble 
constant Hq = 100 h kms™ 1 Mpc -1 , and assuming n(z) and u{z) to be independent of redshift, At 
measured between redshifts 0.35 and 0.20 constrains no < 2 x 10 -4 h Mpc -1 at 95% confidence. 

A naive calculation of At using Equation (|10p for the concordance ACDM model (CIm 



0.258, fi A = 0.742) obtained from the analysis of the 5-year WMAP data (|Dunklev et all 120081 ) 
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yields At = —0.30 ± 0.26 at 95% confidence. This shows that there is a slight tension between the 
results of current measurements of the BAF and of the SNela under the currently accepted world 
model. More generally, a similar tension, i.e., a brightening of the SNe, between measurements 
of the cosmological parameters b y using standard ru l ers and standard c andles has been reported 
before (|Bassett fc Kuna l2004al lbl; IPercival et all 120071 ; lLazkoz et al.l 120081). SNe br ightening is not 
impossible in models that in volve axion-photon mixing (jBassett fc Kuna l2004bl ) or chameleon- 
photon mixing (|Burragell2008l ) if the corresponding particles are abundantly produced during SNela 
explosions. However, a negative value of At could also indicate the presence of a systematic bias 
in the distance measurements based upon the SNela brightness or the BAF, e.g., overcorrection for 
extinction in the host ga laxy of the SNela brightnesses or magnification bias in the SNela selection 
(jWilliams &: Songj 12004! ) . Note that the prior, At > 0, improves the magnitude of the uncertainty 
on At (from 0.26 to 0.13). The 95% confidence interval shrinks with the prior because we sample 
only from the rapidly falling tail of the posterior. 



3. Future Constraints 

In the future, the constraints from both the SNela and the BAF observations will improve 
in accuracy and will cover a wider redshift range. The Baryon Oscillation Spectroscopic Survey 
(BOSS) is currently underway and plans to measure the BAF in luminous red galaxies at redshifts 
z = 0.35 and z = 0.6. The key improvements would be the larger redshift range and the power 
to resolve the BAF both in the line-of-sight direction (constrains H) and the transverse direction 
(constrains -Da). This would remove the weak world-model dependence in our present analysis. 
The angular diameter distances to these redshifts would be measured to an accuracy of ~ 1% 
(http://www.sdss3.org/). In parallel, the Supernovae Legacy Suryey (S NLS), when complete, ex- 



pects ~ 700 SNela in the redshift range < z < 1.7 (jAstier et aZ.1 120061 ) . The uncertainty on the 
estimate of the distance moduli to redshifts z = 0.35 and z = 0.6 will be roughly four times better 
with the increased numbers. Using the test of the duality relation described above, At between 
z = 0.35 and z = 0.6 would be constrained to better than 0.07 (95% confidence), independent of 
the adopted cosmological model. The constraint on nu would become na < 5.4 x 1CP 5 h Mpc -1 . 

BOSS will also use the Lya forest in the spectra of bright quasi-stellar objects to measure the 
BAF at redshift z ~ 2.5 with an accuracy of ~ 1.5 percent. No current or planned SNela surveys 
expect to detect SNela at such a high redshift. However the highest redshift (~ 1.7) measurements 
of -Dl from the SNLS could potentially be used in conjunction with the -Da measurement to get 
a constraint on the transparency of the universe by marginalizing over different world models. 
Interestingly, there have been recent effor ts to calibra te GRBs as standard candles and to extend 



the Hubble diagram to higher redshifts ([Liang et aZ.I l2008l ) . The SNela at low redshift and the 



GRBs at high redshifts can provide a measurement of the difference between the DM between 
redshifts 0.35 and 2.5. We optimistically assume that the difference in the DM to these redshifts 
can be measured with an accuracy of ~ 0.1 similar to the one obtained from the analysis of SNela at 
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z = 0.2 and z = 0.35 in Section [2j These measurements shall then constrain At between redshifts 
2.5 and 0.35 to an accuracy of 0.2 with 95 percent confidence. This translates into an accuracy on 
no- of ~ 1.1 x 10~ 5 h Mpc -1 . 

In the optimistic future, the uncertainty on D^(z) could, in the absence of damaging sys- 
tematics, diminish arbitrarily as the number of observed SNela grows. However, the precision of 
any BAF measurement is limited by sample variance (the number of independent wavelengths of 
a given fluctuation that can fit in the finite survey volume is limited), ev en when the uncertaint y 



caused by incomplete sampling of the density field (shot noise) is negligible (|Seo &: Eisensteinll2007l ) . 
The sample variance error goes down with the square root of the volume of the survey. To cal- 
culate a representative limit, we consider an optimistic all-sky survey covering the redshift range 
2.45 < z < 2.55. Such a survey can be used to determine Dj±(z = 2.5) to a fractional accuracy of 
~ 0.004 (95% confidence). This will ultimately constrain the optical depth to redshift z = 2.5 to 
r < 0.008 and hypothetical absorbers to na < 4 x 10~ 7 h Mpc -1 . 



4. Discussion 

We have advocated and analyzed the expected future performance of a simple Tolman test or 
test of the Etherington relation, that is, that the luminosity distance is larger than the angular 
diameter distance by two powers (1 + z), using SNela to measure the luminosity distance and 
the BAF to measure the angular diameter distance. We have shown that this test will eventually 
provide very precise measurements of the conservation of photon phase-space density. 

We performed the test with the limited data available at the present day. We used only the 
ratio of distances at redshifts of z = 0.20 and 0.35 to remove uncertainties about the overall scale. 
We find consistency with a Lorentz-invariant, transparent universe. Our results are consistent 
with all other measures of transparency to date. This is in part because they are not extremely 
precise. Our Tolman test also assumes that the measurements of the SNela and of the BAF are 
not affected by systematic biases with magnitudes that are a significant fraction of the magnitudes 
of the uncertainties. Our test is limited by the precision of the BAF measurement and the redshift 
range over which it has been measured. As we have shown, experiments planned and underway 
will increase the redshift range and improve the overall precision by an order of magnitude. 

The most precise transparency measurements at visible wavelengths today are statistical angu- 
lar difference measurements, which can only constrain atte nuation correlated w ith specific types of 



absorbing struct ures in the universe (e.g., Mgll absorbers. iMenard et a/.ll2008l : clusters of galaxies, 



Bovy et al.ll20 08). The simple Tolman test performed here limits the full, unclustered, line-of-sight 



attenuation between two redshifts. 

The technique used in this paper provides a test of transparency that is not very sensitive to 
astrophysical assumptions, both because the BAF has a straightforward origin during an epoch 
in which growth of structure is linear and the dominant physics is well understood, and because 
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there is no significant "evolution" with cosmic time for which we must account. This is in contrast 
to other methods for measuring angular diameter distances and brightnesses, where there are 
no precisely "standard" rulers, and evolution is dramatic with redshift. On the other hand, the 
ultimate precision of any test of this type may come from the finite comoving volume in the 
observable universe. Cosmic variance will dominate the BAF error budget eventually. 

The SNela samples have been corrected as best as possible for line-of-sight extinction by fitting 
an empirical correlation of extinction with a change in color. However, t here are a few problems 
with this approach. First, this approach cannot correct for "gray" dust (|Aguirrelll999l ). Second, 
this approach can also not correct for a monopole component; it only corrects for components 
that show variations around the mean level. Third, these corrections will be wrong or fooled if 
there are intrinsic relationships between color and luminosity for SNela. Fourth, the empirical 
corrections found by these projects tend to be odd i n the context of what is expected from the 



reddening and attenuation b y dust in the Milky Way (|Joeveerlll982l : IConley et a/Jl2007l ; Ellis et al 



2008 



Nobili &: Goobarl 120081 ) . The Tolman test is sensitive to any kind of absorber and makes no 



assumptions about the wavelength dependence or fluctuations of the opacity. Given that the SNela 
results have been corrected for a color-brightness relation, the test presented here looks at the 
mean opacity toward SNela of the fiducial color to which the compiled SNela have been corrected. 

The best-fit value of At obtained from our analysis is negative, i.e., SNela are brighter than 
expected from the angular diameter distance measurements using the BAF. A conversion of dark 
sector particles into photons could provide a physical explanation for this result. However, this 
could also indicate the presence of a systematic bias in either the SNela or BAF experiments. A 
test such as the one presented in this paper is a useful tool to bring such biases to light. 

At present, because the differences among competitive world models are not large over the 
redshift range 0.20 < z < 0.35, our test is not yet sensitive enough to rule out extreme axion or 
"gray" dust models that reconcile SNela results with an Einstein-de Sitter universe by using effective 
opacity to adjust the inferred redshift -luminosity-distance relation. Howeve r, these models will all 
be severely constrained within the next few years (see also ICorasanitil 120061 ) . 
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Fig. 1. — Distanc e- modulus-redshif t relation. Filled black squares with uncertainty bars show the 
SNela data (from lDavis et all 120071 ) used in Samples A (left panel) and B (right panel). Open red 
squares show the distance moduli £>M(0.20) = 40.14 ± 0.06 and DM(0.35) = 41.48 ± 0.07 
confidence) inferred from the fits to the data. 
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Fig. 2. — Posterior distribution of At between z = 0.35 and z = 0.20 obtained from the Bayesian 
analysis described in Section [2j The 68, 95 and 99% confidence upper limits are indicated by the 
corresponding dashed lines. 



